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Abstract 

There is a growing interest in the cosmologists for theories with 
negative energy scalar fields and creation, in order to model a re- 
pulsive gravity. The classical steady state cosmology proposed by 
Bondi, Gold and Hoyle in 1948, was the first such theory which used 
a negative kinetic energy creation field to invoke creation of matter. 
We emphasize that creation plays very crucial role in cosmology and 
provides a natural explanation to the various explosive phenomena 
occurring in local (z < 0.1) and extra galactic universe. We exemplify 
this point of view by considering the resurrected version of this theory 
- the quasi-steady state theory, which tries to relate creation events 
directly to the large scale dynamics of the universe and supplies more 
natural explanations of the observed phenomena. 

Although, the theory predicts a decelerating universe at the present 
era, it explains successfully the recent SNe la observations (which 
require an accelerating universe in the standard cosmology), as we 
show in this paper by performing a Bayesian analysis of the data. 
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1. INTRODUCTION 

Remarkable progress has been made in various types of astrophysical and 
cosmological observations in recent years. Among these, the accurate mea- 
surements of the anisotropics in the CMB made by the WMAP experiment 
appear to offer the most promising determination of the cosmological pa- 
rameters. The results of the WMAP experiment are however often quoted as 
providing a direct evidence for an accelerating universe, which is though not 
correct. The cosmological constraints as established by the WMAP team 
(Spergel et al. 2003) entirely rely on the power law spectrum assumption 
and could be erroneous (Kinney 2001; Hannestad 2001). Taken on their 
face value, the WMAP observations are fully consistent with the deceler- 
ating models like the CDM Einstein-de Sitter model (Vishwakarma 2003; 
Blanchard 2005). 

The possibility of an accelerating universe in fact emerges from the mea- 
surements of distant SNe la, which look fainter than they are expected in the 
standard decelerating models. This observed faintness is generally explained 
by invoking some hypothetical source with negative pressure generally known 
as 'dark energy'. This happens because the metric distance of an object, out 
to any redshift, can be increased by incorporating a 'fluid' with negative pres- 
sure in Einstein's equations and hence the object looks fainter. The simplest 
and the most favoured candidate of dark energy is a positive cosmological 
constant A, which is though plagued with the horrible fine tuning problems 
- an issue amply discussed in the literature. This has led a number of cos- 
mologists to resort to scalar field models called quintessence whose function 
is to cause the scale factor to accelerate at late times by violating the strong 
energy condition. While the scalar field models enjoy considerable popular- 
ity, they have not helped us to understand the nature of dark energy at a 
deeper level. By and large, the scalar field potentials used in the literature 
have no natural field theoretical justification and have to be interpreted as 
a low energy effective potential in an ad-hoc manner. Moreover they also 
require fine tuning of the parameters in order to be viable (to find several 
other shortcomings, see for example, Padmanabhan 2005). 

As desperate times call for desperate measures, the cosmologists, in order 
to model the dark energy, have now turned to 'phantom' or 'ghost' scalar 
field models with negative kinetic energy (Caldwell 2002; Carroll, Hoffman & 
Trodden 2003; Gibbons 2003; Singh, Sami & Dadhich 2003; Sami & Toporen- 
sky 2004). The classical steady state cosmology proposed by Bondi, Gold and 
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Hoyle in 1948, was the first such theory which used a negative kinetic energy 
creation field to invoke creation of matter. It is interesting to note that, dis- 
tinct from all the existing big bang models at that time, this model predicted 
an accelerating universe. However, it is unfortunate that the theory was not 
given any credit (which it deserved, despite the difficulties associated with 
it) when the SNe la observations started claiming an accelerating universe 
in 1998. 

Once the cosmologists have lost their inhibitions about negative energy 
fields, the time is ripe for considering the idea that the creation of matter 
plays important role in cosmology. We exemplify this point of view by con- 
sidering the resurrected version of the classical steady state theory, namely 
the quasi-steady state cosmology (QSSC) which has not been given proper 
attention as it deserves. This theory was proposed by Hoyle, Burbidge and 
Narlikar in 1993 (1993; 1995), wherein the introduction of negative kinetic 
energy scalar field is not ad-hoc but is required to ensure that matter creation 
does not violate the law of conservation of matter and energy. However, first 
we emphasize that the idea of creation of matter is already present in general 
relativity, though hidden behind some simplifying assumptions. 

With a suitable Lagrangian for the source terms, the Einstein field equa- 
tions can be written as 



Rik — ^dikR — —8nG 



^(matter) rp(A) rp(4>) . 
1 ik ' ik ' ik 



where we have considered the speed of light c = 1. The only constraint 
on the source terms, which is imposed by this equation, is the conservation 
of right hand side through the Bianchi identities: [Rij — ^Rg^'i = = 
^(matter) rp(A) rp{4>) _|_ ^ J ^ i m piyi n g that only the sum of all the energy- 
momentum tensors is conserved, individually they are not. If we take them 
conserved separately, as is in practice among the cosmologists, it can be done 
only through the additional assumption of no interaction (minimal coupling) 
between different source fields, which though seems ad-hoc and nothing more 
than a simplifying assumption. On the contrary, interaction is more natural 
and is a fundamental principle. Of course some ideal cases are consistent 
with the idea of minimal coupling, for example, j^™ attcr ) with a constant 
A. However, imposing this assumption on non-trivial cases would result in 
losing some important information. For example, taken on the face value, 
a time-dependent A, with matter, implies matter creation and results in a 



3 



Machian model (Vishwakarma 2002a). However, if one makes an additional 
assumption of no interaction between 7;(™ atter ) an( j A(t), these features are 
lost and A(t) reduces to a constant. It is well known that even if one considers 
the Robertson- Walker spacetime (to avoid non-Machian Godel's solution of 
Einstein field equations), there still exists a non-Machian solution of Einstein 
field equations - the de Sitter solution. Creation has many more attractive 
features. It has been shown how the scalar creation field helps in resolving 
the problems of singularity, flatness and horizon in cosmology (Narlikar & 
Padmanabhan 1985). Such a negative energy creation field is responsible for 
a non-singular bounce from a high non-singular density state, as has been 
shown by Hoyle and Narlikar (1964). This idea has been recently used by 
Steinhardt and Turok in their oscillatory model (Steinhardt & Turok 2002). 
The quasi-steady state cosmology (QSSC) is also a Machian theory which 
is derived from an action principle based on Mach's Principle, and assumes 
that the inertia of matter owes its origin to other matter in the universe. 
The stress-energy tensor for creation (corresponding to T^' in equation (1)) 
is given by 

^creation = _ f + 1^ j (2) 

where / is a positive coupling constant and the gradient Cj = d(f>/dx l is the 
contribution from a trace-free zero rest mass scalar field <p of negative energy 
and stresses. The A in this theory (corresponding to a = —hg i j t /%'KG 
of (1)) appears as a constant of nature with its value ~ — 2 x 10~ 56 cm" 2 , 
which falls within the normally expected region of the magnitude of the 
cosmo logical constant. Note however that its sign is negative, which is a 
consequence of the Machian origin of the cosmological constant. The theory 
does not face the cosmological constant problem mentioned earlier. In fact, 
the A in the QSSC does not represent the energy density of the quantum 
fields, as this model does not experience the energy scales of quantum gravity 
except within the local centres of creation. The theory offers a purely stellar- 
based interpretation of all observed nuclei including the light ones (Burbidge 
et al. 1957; Burbidge & Hoyle 1998). In the following, we demonstrate in 
brief the main features of this cosmology and how it confronts the various 
observations (for more details, one can consult (Sachs, Narlikar & Hoyle 1996; 
Hoyle, Burbidge & Narlikar 2000)). 

The QSSC represents a cyclic universe with its Robertson- Walker scale 



4 



factor given by 

S(t) = e*/ p [l + 77 cos(27rr/Q)], (3) 

where the time scales P ~ 10 3 Gyr ^> Q ^ 40 — 50 Gyr are considerably 
greater than the Hubble time scale of 10 — 15 Gyr of the standard cosmology. 
The function r(t) is very nearly like the cosmic time t, with significantly 
different behaviour for short duration near the minima of the function S(t). 
The parameter rj has modulus less than unity, thus preventing the scale factor 
from reaching zero. Typically, r) ~ 0.8 — 0.9. Hence there is no spacetime 
singularity, nor a violation of the law of conservation of matter and energy, 
as happens at the big bang epoch in the standard cosmology. The model has 
cycles of expansion and contraction (regulated respectively by the creation 
field and the negative A) of comparatively shorter period (Q) superposed 
on a long term (P) steady state-like expansion. Creation of matter, which 
occurs through explosive processes, is also periodic, being confined to pockets 
of strong gravitational fields around compact massive objects and the nuclei 
of existing galaxies. Such processes take place whenever the energy of the 
creation field quantum rises above a threshold energy, which is equal to the 
restmass energy of the created Planck particle. 

The model provides a natural explanation to the various explosive phe- 
nomena occurring in local (z < 0.1) and extra galactic universe. By the early 
1960s it had become clear that very large energy outbursts are taking place 
in the nuclei of galaxies. In the decades since then it has been found that 
many active nuclei are giving rise to x-rays, and relativistic jets, detected in 
the most detail as high frequency radio waves. A very large fraction of all of 
the energy which is detected in the compact sources is non-thermal in origin, 
and is likely to be incoherent synchrotron radiation or Compton radiation. 
In addition to this we see several other explosive phenomena in the universe, 
such as jets from radio sources, gamma ray bursts, X-ray bursters, QSOs, 
etc. Generally it is assumed that a black hole plays the lead role in such an 
event by somehow converting a fraction of its huge gravitational energy into 
large kinetic energy of the 'burst' kind. In actuality however, we do not see 
infalling matter that is the signature of a black hole. Rather we see outgoing 
matter and radiation, which agrees very well with the idea of creation events 
formulated in the framework of the QSSC. 

There are several free parameters in the model which are estimated from 
the observations and provide a decelerating universe at the present cycle of 
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expansion. It is then interesting to see how the model explains the SNe la 
and other observations! This is shown in the following. 

2. The High Redshift Supernovae la 

It is generally accepted that metallic vapours are ejected from the SNe explo- 
sions which are subsequently pushed out of the galaxy through pressure of 
shock waves (Hoyle & Wickramasinghe 1988). Experiments have shown that 
metallic vapours on cooling, condense into elongated whiskers of ps 0.5 — 1 mm 
length and psIO -6 cm cross-sectional radius (Donn & Sears 1963; Nabarro & 
Jackson 1958). It can be shown that the extinction from the whisker dust 
adds an extra magnitude Sm(z) to the apparent magnitude m(z) (arising 
from the cosmological evolution) of the SN light emitted at the epoch of 
redshift z, which is given by 

fZ Ay 1 

5m{z) = 1.0857 x k p g0 / (l + z 'f——. (4) 

Jo H\ z ) 

where k is the mass absorption coefficient and p go is the whisker grain density 
at the present epoch. The net apparent magnitude is then given by 

m QCt (z) = m(z) +5m(z). (5) 

By taking account of this effect, it has been shown that this kind of dust 
extinguishes radiation travelling over long distances and decelerating models 
without any dark energy (for example, the Einstein-de Sitter model) can also 
explain high redshift SNe la observations successfully (Vishwakarma 2002b; 
2003; 2005). QSSC in fact resorts to this dust to explain not only the SNe 
la observations but also CMB, as we shall see in the following. It has been 
shown (Narlikar, Vishwakarma & Burbidge 2002; Vishwakarma & Narlikar 
2005) that by taking account of this effect, QSSC explains successfully the 
SNe la data from Perlmutter et al. (1999) and also shows an acceptable 
fit to the 'gold sample' of 157 SNe la recently published by Riess et al. 
(2004) which, in addition to having previously observed SNe, also includes 
some newly discovered highest-redshift SNe la by the Hubble Space Tele- 
scope. Though this sample is believed to have a 'high-confidence' quality of 
the spectroscopic and photometric record for individual supernovae, we note 
that there are some SNe (1997as, 1997bj, 2000eg, 2001iw, 2001iv) in this 
sample which do not seem to be consistent with any of the models generally 
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considered in the fitting and appear as general outliers (see the encircled 
SNe in Fig. 1 and Fig. 2). By excluding these points, the fit to different 
models improves considerably. For example, the x 2 value per degrees of free- 
dom (dof) for the best-fitting QSSC model reduces to 1.18 from the earlier 
X 2 /dof=1.30 obtained from the full sample of 157 points (Vishwakarma & 
Narlikar 2005). The fit to the standard (flat ACDM) cosmology improves 
tremendously from x 2 /dof=1.14 (from 157 points) to % 2 /dof=0.99 (from 152 
points). The details of the fit (in the case of the frequentist approach) can 
be found in (Vishwakarma & Narlikar 2005). 
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Figure 1: Some best-fitting models are compared with the 'gold sample' of SNe 
la data with 157 points as considered by Riess et al. (2004). The solid curve cor- 
responds to the QSSC model with the whisker dust, the dotted curve corresponds 
to the flat ACDM model, the dashed curve corresponds to the spherical ACDM 
model, and the dashed-dotted curve corresponds to the Einstein-de Sitter model. 
The models differ significantly for z > 1.2. The encircled points seem to be general 
outliers which are missed by all the models. 

Though there is not a clearly denned value of x 2 /dof for an acceptable fit, 
a 'rule of thumb ' for a moderately good fit is that x 2 should be roughly equal 
to the number of dof. A more quantitative measure for the goodness- of- fit is 
given by the ■^-•probability. If the fitted model provides a typical value of x 2 
as x at n dof, this probability is given by 



Roughly speaking, it measures the probability that the model does describe 
the data and any discrepancies are mere fluctuations which could have arisen 
by chance. To be more precise, Q(x,n) gives the probability that a model 
which does fit the data at n dof, would give a value of x 2 as large or larger than 
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Figure 2: Modified Hubble diagram of the 'gold sample' of SNe la minus a fiducial 
model (f2 m o = 0, Qao = 0). The relative magnitude(Am net = m net — mgduciai) is 
plotted for some best-fitting models, by using the original error bars. The solid 
curve corresponds to the QSSC model with the whisker dust, the dotted curve 
corresponds to the flat ACDM model, the dashed curve corresponds to the spherical 
ACDM model, and the dashed-dotted curve corresponds to the Einstein-de Sitter 
model. The encircled points seem to be general outliers which are missed by all 
the models. 



x. If Q is very small, the apparent discrepancies are unlikely to be chance 
fluctuations and the model is ruled out. It may however be noted that the x 2 - 
probability strictly holds only when the models are linear in their parameters 
and the measurement errors are normally distributed. It is though common, 
and usually not too wrong, to assume that the x 2 ~ distribution holds even 
for models which are not strictly linear in their parameters, and for this 
reason, the models with a probability as low as Q > 0.001 are usually deemed 
acceptable (Press et al. 1986). Models with vastly smaller values of Q, say, 
10~ 18 are rejected. The probability Q for the best-fitting QSSC to the full 
sample is obtained as 0.007, which is though very small, but acceptable. 
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By excluding the above-mentioned 5 outliers, Q improves to 0.062. The 
corresponding probabilities in the case of the standard ACDM cosmology 
are obtained as 0.109 and 0.534. 

We note that the fit to the QSSC is considerably worse than those in the 
standard ACDM cosmology. However, one cannot compare the relative mer- 
its of the models on the basis of the x 2 -probability (frequentist approach), 
which uses the best-fitting parameter values and hence judges only the max- 
imum likely performance of the models. The more appropriate theory for 
such comparisons is the Bayesian theory which does not hinge upon the 
best-fitting parameter values and evaluates the overall performance of the 
models by using average likelihoods (rather than the maximum likelihoods), 
given by the Bayes factor B. The theory employs the premise that if we 
assume an equal prior probability for competing models, the probability for 
a given model is proportional to the marginalised likelihood called evidence. 
We have described this theory, in brief, in the Appendix (for more details, 
see (Drell, Loredo & Wasserman 2000; John & Narlikar 2002)). 

In order to calculate the Bayes factor B for the two models QSSC and 
the standard ACDM, first we have to fix the prior probabilities for the free 
parameters. While the flat QSSC has four free parameters kp^Hq 1 , Qao, 
z max and M. (see the Appendix of (Vishwakarma & Narlikar 2005)), the 
standard ACDM has only two free parameters Sl m0 and Ai. We would like 
to mention that the whisker dust was already introduced in the QSSC in 
order to explain the CMB which put a constraint on the density of the dust 
p g o ~ 10~ 34 g cm~ 3 (Narlikar et al. 2003). This value, taken together with 
the observational constraints on k (Wickramasinghe & Wallis 1996) and H 
(Freedman & Turner 2003) from other observations, supply a value of the 
parameter np^H^ 1 close to the best-fitting value estimated from the SNe la 
observations. Hence we assign a prior on the parameter npgoH^ 1 that it lies 
in the range np^H^ 1 G [3.5,6]. We also note that ^ao m the QSSC does 
not receive any significant contribution from £7ao < —0.3. Taking account of 
this and the theoretical constraint that A in the QSSC is negative, we assign 
a prior on the parameter f2Ao G [—0.3, 0]. To the rest two parameters in the 
QSSC, about which we do not have prior information, we assign liberal priors: 
<2max £ [5, 10] (to be consistent with the highest redshift « 7 observed so far) 
and Ai G [41,45] (which is the common parameter). For the parameter Q m0 
in the flat ACDM model, we assume that Sl m0 G [0, 1] (which is equivalent to 
assigning f^Ao £ [0, !])• It may be noted that the fl\o in the two models are 
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altogether different quantities (though they have been denoted by the same 
symbol in order to match the general convention) and there is no reason to 
assign the same probability for them in the two different models. 

When calculated for the full 'gold sample' of 157 points, these prior prob- 
abilities give a Bayes factor favouring the standard ACDM over the QSSC as 
B = 2.77, which though indicates an evidence against the QSSC, however, 
the evidence is not definite and is not worth more than a bare mention (for 
the interpretation of B, see the Appendix). Our assigning f^o G [—0.3, 0] in 
the QSSC can raise eyebrows, as this probability is very conservative com- 
pared with the one in the ACDM. However, assigning this probability is due 
to the reason that the likelihood for A in the QSSC does not receive any 
significant contribution from Q\o < —0.3, as mentioned earlier. For exam- 
ple, increasing the domain of A in its prior to Qao G [—1,0] in the QSSC, 
results in lowering the likelihood of the model, as expected. This gives a 
Bayes factor B = 9.30, which indicates that the evidence against the QSSC 
is definite, though not strong. 

One should also note that a proper assessment of the probability for a 
model is given by p = 1/(1 + 5). Thus for the above-mentioned two choices 
of the prior probabilities, the corresponding probabilities for the QSSC are 
0.27 and 0.10, which are reasonably good probabilities. 

It is interesting to note that the whisker dust, which is a vital ingredient of 
the QSSC, does not make any significant improvement in the fit to the ACDM 
cosmology. It may be argued that this kind of dust can create too much 
optical depth for the high redshift objects and they need to be excessively 
bright in order to be seen. However, from our calculations, we find that 
the objects over the present cycle right upto the maximum redshift will be 
fainter, at the most, by ~6 magnitudes only and it is possible to see even 
sources from many previous cycles, though they will be very faint. 

In Fig. 1, we have compared some best-fitting models with the actual 
data points. In order to have a better visual comparison of these models, 
we magnify their differences by plotting the relative magnitude with respect 
to a fiducial model f2 m0 = ^ao = 0, without any whiskers (which has a 
reasonably good fit: x 2 /dof = 1.2). This has been shown in the 'modified' 
Hubble diagram in Fig. 2. In the following we describe, in brief, the other 
important features of the QSSC. 

3. The CMB 
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As far as the origin and nature of the CMB is concerned, the QSSC uses a fact 
that is always ignored by standard cosmologists. If we suppose that most of 
the 4 He found in our own and external galaxies (about 24% of the hydrogen 
by mass) was synthesized by hydrogen burning in stars, the energy released 
amounts to about 4.37 x 10~ 13 erg cm~ 3 . This is almost exactly equal to 
the energy density of the microwave background radiation with T = 2.74°K! 
In the standard cosmology, this has to be dismissed as a coincidence. Thus 
according to the QSSC, the CMB is the relic starlight left by the stars of 
the previous cycles which has been thermalized by the metallic whisker dust 
emitted by the supernovae. As a typical cycle proceeds from the maximum of 
scale factor towards the next minimum, the wavelengths of the starlight from 
the previous cycle are shortened since the universe contracts by a considerable 
factor. It has been shown (Narlikar et al. 1997) that at wavelengths 100 m-20 
cm, sufficient optical depth exists for this radiation to thermalize in about 
twenty cycles. The production of microwaves in this fashion goes on in each 
cycle, and the process of frequent absorption and re-radiation by whiskers will 
eventually generate a uniform background, except for the contribution from 
the latest generation of clusters. These will stand out as inhomogeneities on 
the overall uniform background arising from certain intrinsic inhomogeneities 
of the process as well as from the cosmological model. A quantitative analysis 
shows that this process requires an intergalactic dust of density of ~ 10 -34 
g cm -3 , which is very close to the best-fitting value estimated from the SNe 
la observations. The theory also explains the peaks at I ~ 200 and / ~ 600 
which are related, in this cosmology, to the clusters and groups of clusters 
(Narlikar et al. 2003). Also, we have taken stock of the WMAP observations 
in (Narlikar, Burbidge & Vishwakarma 2007). 

Though these studies do not give predictions as sharp as those given by 
the standard big bang cosmology, however, one should note the attitudinal 
difference between this approach and the standard one. In the big bang cos- 
mology, the inferences are related to the postulated initial conditions prevail- 
ing well beyond the range of direct observations (at redshifts 1100). Whereas 
the QSSC interpretation links the inhomogeneities of the radiation field to 
those of the matter field, on which we do not have very accurate data at 
present, but which may be observable one day. 

It may also be worthwhile to mention that there are claims that like the 
dipole, the quadrupole and the octopole harmonics of the CMB spectrum 
also have their origin in the solar system (Starkman et al. 2004). If this is 
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correct, then subtracting this foreground contribution from the rest of the 
signal (in order to have the temperature fluctuations only at the time of the 
big bang) would render the inflationary model in serious trouble. 

4. The Non-Baryonic Dark Matter 

Unlike the standard big bang cosmology, the QSSC allows the dark matter 
to be baryonic. It may be recalled that the standard cosmology predicts the 
existence of non-baryonic, though as yet undetected, particles to solve the 
problems of structure formation and of the missing mass in bound gravita- 
tional systems such as galaxies and clusters of galaxies. The most favoured 
candidate of non-baryonic dark matter postulated by many astrophysicists, 
cosmologists and particle physicists is a massive but very weakly interacting 
particle called WIMP (Weakly Interacting Massive Particle), a hypothetical 
elementary particle that was produced moments after the Big Bang. Cur- 
rently there are a number of WIMP detection experiments underway. Among 
these, the DAM A experiment (Bernabei et al. 2003), which measures the 
annual modulation in WIMP interactions with the sodium-iodide detectors 
caused by the earth's rotation around the sun, is the only one to have claimed 
a positive signal. However, the results of this experiment are controversial as 
other more sensitive searches have not detected nuclear recoils due to WIMP 
interactions (Akerib et al. 2004; Angloher et al. 2005) and concluded that al- 
most all the events measured by DAMA were from neutrons, and should not 
be attributed to scattering events from dark-matter WIMPs. It is therefore 
fair to say that this scheme has still to demonstrate its viability. 

However, in the framework of the QSSC, the dark matter need not be 
necessarily non-baryonic. It can be in the form of baryonic matter being the 
relic of very old stars of the previous cycles. A typical QSSC cycle has a 
lifetime long enough for most stars of masses exceeding ~ 0.5 — O.7M to 
have burnt out. Thus stars from previous cycles will be mostly extinct as 
radiators of energy. Their masses will continue, however, to exert a gravi- 
tational influence on visible matter. The so-called dark matter seen in the 
outer reaches of galaxies and within clusters may very well be made up, at 
least in part, of these stellar remnants. 

It may be timely to mention that the recent data on distant x-ray clusters 
obtained from XMM and Chandra projects indicate that the observed abun- 
dances of clusters at high redshift, taken at face value, give 0.9 < Q m0 < 1.07 
(at 1 a) (Blanchard 2005). This favours a matter-dominated model and is 
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consistent with the value of f2 m o in the QSSC estimated from the SNe la and 
CMB observations. However, it is hard to reconcile with the concordance 
model. 

5. Conclusion 

In order to explain the current observations in the framework of the stan- 
dard cosmology, one has to trust upon a preposterous composition for the 
constituent of the universe which defies any simple explanation, thereby pos- 
ing probably the greatest challenge theoretical physics has ever faced. We 
think this is the right time to consider seriously alternative theories which 
present more natural explanations to the observed phenomena, especially 
when there is no independent observational evidence for non-baryonic dark 
matter, dark energy and inflation, neither they have a firm basis in a well 
established theory of particle physics. Furthermore it is always necessary for 
healthy science to have an alternative model to the dominant paradigm. 

Acknowledgement RGV thanks the Abdus Salam ICTP, for providing hos- 
pitality and travel assistance (under his associateship programme), where a 
part of the work was done. 

APPENDIX 

The Bayes factor is a ratio of average likelihoods (rather than the maximum 
likelihoods used for model comparison in frequentist statistics) for two models 
Mj, and Mj, and is given by 

C(M t ) _ p{D\M % ) 
~~ C(Mj) ~ p(D\MjY {AA) 

where the likelihood for the model Mj, £(Mj) is the probability p(D\Mi) to 
obtain the data D if the model Mj is the true one. For a model Mj with free 
parameter, say, a and (5 (generalization for the models with more parameters 
is straight forward), this probability is given by 

C{Mi) =p(D\Mi) =JdaJ dp P (a\Mi)p(f3\Mi)jCi(a,P), (A.2) 

where p(a|Mj) and p{(5\Mi) are the prior probabilities for the parameters 
a and (5 respectively, assuming that the model Mj is true. £j(a,/9) is the 
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likelihood for a and (3 in the model Mj and is given by the usual x 2 -statistic: 



£i(a,f3) = exp 



(A.3) 



For a flat prior probabilities for the parameters a and ft, i.e., assuming that 
we have no prior information regarding a and (5 except that they lie in 
some range [a, a + Act] and \J3, (3 + A/3], we have p(a|Mj) = 1/Act and 
p(f3\Mi) = l/A/3. Hence the expression for the likelihood of the model Mj 
reduces to 



C(Mi) = 



1 1 

Aa A(3 J a 



ra+Aa rf3+A/3 




/ / exp 




la J/3 


2 



dp da. 



(A.4) 



The Bayes factor Bij, given by (A.l), which measures the relative merits of 
model Mi over model Mj, is interpreted as follows (Drell, Loredo & Wasser- 
man 2000; John & Narlikar 2002; and the references therein). If 1 < B^ < 3, 
there is an evidence against Mj when compared with Mj, but it is not worth 
more than a bare mention. If 3 < B^ < 20, the evidence against Mj is 
definite but not strong. For 20 < B^ < 150, this evidence is strong and for 
B^ > 150, it is very strong. 
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